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Abstract We use an early twentieth century (1908-1958) 
atmospheric reanalysis, based on assimilation of surface 
and sea level pressure observations, to contrast atmospheric 
circulation during two periods of persistent drought in 
North America: 1932-1939 (the ‘Dust Bowl’) and 1948- 
1957. Primary forcing for both droughts is believed to 
come from anomalous sea surface temperatures (SSTs): a 
warm Atlantic and a cool eastern tropical Pacific. For 
boreal winter (October-March) in the 1950s, a stationary 
wave pattern originating from the tropical Pacific is pres- 
ent, with positive centers over the north Pacific and north 
Atlantic ocean basins and a negative center positioned over 
northwest North America and the tropical/subtropical 
Pacific. This wave train is largely absent for the 1930s 
drought; boreal winter height anomalies are organized 
much more zonally, with positive heights extending across 
northern North America. For boreal summer (April-Sep- 
tember) during the 1930s, a strong upper level ridge is 
centered over the Great Plains; this feature is absent during 
the 1950s and appears to be linked to a weakening of the 
Great Plains low-level jet (GPLLJ). Subsidence anomalies 
are co-located over the centers of each drought: in the 
central Great Plains for the 1930s and in a band extending 
from the southwest to the southeastern United States for the 
1950s. The location and intensity of this subsidence during 
the 1948-1957 drought is a typical response to a cold 
eastern tropical Pacific, but for 1932-1939 deviates in 
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terms of the expected intensity, location, and spatial extent. 
Overall, circulation anomalies during the 1950s drought 
appear consistent with the expected response to the 
observed SST forcing. This is not the case for the 1930s, 
implying some other causal factor may be needed to 
explain the Dust Bowl drought anomalies. In addition to 
SST forcing, the 1930s were also characterized by massive 
alterations to the land surface, including regional-scale 
devegetation from crop failures and intensive wind erosion 
and dust storms. Incorporation of these land surface factors 
into a general circulation model greatly improves the 
simulation of precipitation and subsidence anomalies dur- 
ing this drought, relative to simulations with SST forcing 
alone. Even with additional forcing from the land surface, 
however, the model still has difficulty reproducing some of 
the other circulation anomalies, including weakening of the 
GPLLJ and strengthening of the upper level ridge during 
AMJJAS. This may be due to either weaknesses in the 
model or uncertainties in the boundary condition estimates. 
Still, analysis of the circulation anomalies supports the 
conclusion of an earlier paper (Cook et al. in Proc Natl 
Acad Sci 106:4997, 2009), demonstrating that land deg- 
radation factors are consistent with the anomalous nature of 
the Dust Bowl drought. 
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1 Introduction 

Recent years have seen rapid advances in our understand- 
ing of the dynamical mechanisms governing global 
hydroclimatic variability (Hoerling and Kumar 2003; 
Schubert et al. 2009; Seager et al. 2003, 2005a; Seager 
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2007). Modeling and empirical studies have revealed 
robust teleconnections between variations in sea surface 
temperatures (SSTs) and precipitation and temperature 
patterns in distant regions, highlighting the possibility for 
improved prediction of climate on decadal timescales 
(Keenlyside et al. 2008). North America (NA) is one 
region where there has been a great degree of success 
simulating hydroclimatic variability using SST forcing. 
Much of the variance and persistence in drought and plu- 
vial events over NA can be dynamically linked to SST 
anomalies and their associated climate patterns originating 
in the tropical Pacific (El Nino Southern Oscillation, 
ENSO) and the north Atlantic (Atlantic Multidecadal 
Oscillation, AMO). 

On interannual to decadal timescales, ENSO is the 
major global driver of hydroclimatic variability in the 
tropics, subtropics, and midlatitudes, including NA 
(Herweijer et al. 2006, 2007; Seager et al. 2005b; Seager 
2007). When the eastern tropical Pacific is colder than 
normal (i.e., La Nina phase of ENSO), precipitation over 
southern NA is suppressed, with the drying arising from 
from two separate, but ultimately complementary, 
dynamical mechanisms. The first is a zonally and hemi- 
spherically symmetric mean meridional circulation 
response, driven by pan-tropical cooling originating from 
the cold tropical Pacific SSTs (Seager et al. 2003, 2005a). 
Cooler tropics lead to a reduction in equator- to-pole tem- 
perature and pressure gradients and a poleward contraction 
of the westerly jets. Transient eddies penetrate less deeply 
into the tropics and deposit their momentum in the mid- 
latitudes, forcing mass convergence in the upper atmo- 
sphere which, to satisfy continuity, drives anomalous 
descent, suppressing precipitation in midlatitude zonal 
bands in both hemispheres. The second mechanism is a 
stationary wave emanating from the tropical Pacific, often 
associated with the Pacific North America pattern (PNA) 
(Trenberth et al. 1998; Wallace and Gutzler 1981). During 
La Nina, this stationary wave is characterized by positive 
height anomalies in the north Pacific and along the eastern 
coast of NA, with negative height anomalies over the 
subtropical Pacific and extending southeast from Alaska 
and the Pacific Northwest (Trenberth et al. 1998). Both 
mechanisms lead to drying in the southwestern and 
southeastern regions of NA. 

The direct influence of ENSO on NA hydroclimate is 
strongest during boreal winter and is weaker during the 
warm season, but can be enhanced indirectly through soil 
moisture anomalies and Atlantic SSTs. Soil moisture acts 
as a seasonal bridge, providing a source of system memory 
at the land surface and allowing cold season moisture 
anomalies to carry over into the warm season, prolonging 
and often intensifying drought anomalies (Schubert et al. 
2004a; Wang et al. (2010)). The Pacific can also act 


indirectly during the warm season through it’ s influence on 
Atlantic SSTs (Enfield and Mayer 1997), which have a 
stronger impact, relative to ENSO, during the warm season 
(Kushnir et al. 2010; Wang et al. 2010). There is some 
speculation that forcing from the extratropical north Pacific 
may also be an important driver of NA hydroclimatic 
variability, especially on interdecadal time scales, mani- 
festing as the Pacific Decadal Oscillation (PDO) (Mantua 
et al. 1997; Mantua and Hare 2002; McCabe et al. 2004, 
2008). However, the underlying dynamics of Pacific dec- 
adal variability are poorly understood (McCabe et al. 
2008), and the associated atmospheric anomalies look quite 
similar to responses due to ENSO forcing. This has led to 
some speculation that Pacific decadal variability may be 
partially or wholly an extratropical low frequency expres- 
sion of ENSO (Mantua and Hare 2002; McCabe et al. 
2008). 

More recently, evidence has emerged for a significant 
influence of north Atlantic SST variability on NA climate 
(Kushnir et al. 2010; McCabe et al. 2004, 2008; Wang 
et al. 2006). At multidecadal frequencies, this SST vari- 
ability is linked with the AMO which is thought to be 
driven by variations in the strength of the ocean meridional 
overturning circulation (Enfield et al. 2001; Sutton and 
Hodson 2005). When North Atlantic SSTs are warmer than 
normal, precipitation is typically reduced over NA, 
although the specific mechanisms remain elusive (McCabe 
et al. 2004, 2008). As with the north Pacific, SSTs in the 
Atlantic are known to be influenced by ENSO (Enfield and 
Mayer 1997), and so Atlantic variability may not be 
completely independent from ENSO-forced variability 
(Enfield et al. 2001; Sutton and Hodson 2005). The influ- 
ence of north Atlantic SSTs on NA hydroclimate is most 
completely understood during the boreal summer (Enfield 
et al. 2001; Sutton and Hodson 2005, 2007), though recent 
work suggests the Atlantic impact persists throughout the 
year (Kushnir et al. 2010; Wang et al. 2010). 

When forced with observed SSTs, atmospheric general 
circulation models (GCMs) reliably reproduce the inten- 
sity, duration, and spatial expression of most instrumental 
era NA droughts. This includes the Civil War drought 
(1856-1865), the 1870s drought (1870-1877), the 1890s 
drought (1890-1896), the 1950s drought (1948-1957), and 
the most recent turn of the twenty-first century drought 
(1998-2004) (Herweijer et al. 2006; Seager et al. 2005b; 
Seager 2007). The notable exception, however, is the ‘Dust 
Bowl’ drought (1932-1939). While there has been some 
success simulating the onset of the Dust Bowl and the 
general drought pattern (Schubert et al. 2004a, b; Seager 
et al. 2008), few models are able to reproduce the proper 
intensity or location of this drought using SST forcing 
alone. The difficulty may be due to the exceptional land 
degradation that occurred during the Dust Bowl, conditions 
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unique among the instrumental droughts (Hansen and 
Libecap 2004). The land degradation took the form of 
extensive crop failures throughout the Great Plains, leading 
to massive wind erosion and dust storm activity that was 
unprecedented for the historical period. Recent modeling 
work suggests that the feedbacks from the dust aerosols 
and the loss of vegetation played a critical role in driving 
the anomalous temperature and precipitation patterns 
(Cook et al. 2008, 2009; Koven 2006). 

If SST forcing alone is insufficient, and other mecha- 
nisms (such as changes at the land surface) are required to 
explain the anomalous nature of the 1930s Dust Bowl, this 
should be apparent in the dynamical fields, and it will be 
expected that circulation during the 1930s will be notice- 
ably different from circulation during droughts dominated 
by SST forcing, such as the 1950s drought. Up until 
recently, a full dynamical investigation of the 1930s has 
been difficult because of the paucity of atmospheric data 
(other than precipitation and temperature) from this time 
period, although some preliminary analyses have been 
conducted (Bronnimann et al. 2009). Recently, however, a 
new reanalysis product has become available, covering the 
time period from 1908-1958, a span encompassing both 
the 1930s drought (1932-1939) and the 1950s drought 
(1948-1957). For this study, we will use variables from this 
reanalysis (geopotential height, pressure vertical velocity, 
winds, specific humidity) to investigate two questions: (1) 
How different is the atmospheric circulation between the 
1930s (1932-1939) and 1950s (1948-1957) North Ameri- 
can droughts? and (2) How well can the 1930s drought be 
simulated by a GCM forced only by observed SSTs, and is 
land surface degradation required to reproduce the circu- 
lation anomalies and drought pattern? 

2 Data and methods 

In order to facilitate comparisons between the model 
results and the COMPO reanalysis, we chose a common 
baseline period (1921-1929) for the calculation of anoma- 
lies. This period was relatively free of any major persis- 
tent drought or pluvial events, representing a largely 
‘normal’ period in NA hydroclimate, and is independent 
from both drought periods we are interested in investigat- 
ing. It is also a common baseline period available from 
both the model runs (1921-1939) and the reanalysis (1908- 
1958). 

2.1 Early twentieth century reanalysis (1908-1958) 

For this study, we use a recently released atmospheric 
reanalysis for the first half of the twentieth century, 1908- 
1958 (Compo et al. 2006, hereafter, COMPO), a product 


conceptually similar to other data assimilation products, 
such as the NCEP/NCAR reanalysis (Kalnay et al. 1996). 
The COMPO reanalysis is generated using an ensemble 
filter approach. Ensembles of short-term forecasts are 
generated using the state of the art atmospheric general 
circulation model from the NCEP operational Climate 
Forecast System (Saha et al. 2006), operating on an 
irregular Gaussian grid corresponding to a T62 horizontal 
resolution (grid spacing approximately 2° x 2°), with 28 
levels in the vertical, extending to a model top of 0.2 hPa. 
The atmosphere model is driven by time-evolving sea 
surface temperatures and sea ice distributions from the 
HadISSTl.1 dataset (Rayner et al. 2003), which we also 
use to compare and contrast SST patterns during the two 
droughts. The model assimilates observations of surface 
pressure and sea level pressure every six hours from the 
International Surface Pressure Databank version 1.1 and 
ICOADS version 2.4. Variables available from the 
reanalysis at the various pressure levels are geopotential 
heights, temperature, horizontal winds, pressure vertical 
velocities, and specific and relative humidity. 

For our region of interest (NA), the COMPO reanalysis 
compares favorably against independent data and analyses, 
including an independent statistical reconstruction of early 
twentieth century geopotential heights (Griesser et al. 
2008) and gridded precipitation anomalies (Mitchell and 
Jones 2005). The Griesser et al. (2008) product is a sta- 
tistical reconstruction of geopotential height and tempera- 
ture fields from historical upper air and surface data. For 
their reconstruction of geopotential heights, Griesser et al. 
(2008) show high skill over their validation period during 
the winter months, but much less skill during the summer. 
For the full overlapping period and our region of interest 
(1921-1957, averaged over 180°W-40°W and 20°N- 
60°N), geopotential heights between COMPO and Griesser 
et al. (2008) are highly correlated, especially during the 
winter. For October-March (ONDJFM), at 850/500/200 
hPa, the Pearson correlations between the two products are 
0.88/0.80/0.45 (all significant at p < 0.05); for these same 
levels during April-September (AMJJAS) the correlations 
are 0.48/0.38/0.56 (all significant at p < 0.05). Again, we 
note that the Griesser et al. (2008) product has significantly 
lower skill during the warm season, so a reduction in the 
correlation of the two products is not especially surprising. 
Both products also show similar geopotential height 
anomalies during our time periods of interest, including the 
wave train during ONDJFM, 1948-1957 and the upper 
level ridge during AMJJAS, 1932-1939; these features will 
be discussed later. Pressure vertical velocity anomalies 
from COMPO match well with observed precipitation 
anomalies from the CRU 2.1 climate grids (Mitchell and 
Jones 2005) for the two drought periods, with subsidence 
anomalies coincident with negative precipitation anomalies 
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(discussed later). Finally, we note that the COMPO 
reanalysis has already been used, successfully, to look at 
other early twentieth century circulation features (Wood 
and Overland 2009). 

For this study we will be focusing on the NA region, a 
relatively data rich area that should provide more reliable 
estimates than other, data poor areas, such as the Southern 
Ocean or Africa. We use monthly average values of the 
major dynamical variables (geopotential heights, horizontal 
winds, and pressure vertical velocity) and organize the 
analysis into two seasons: boreal autumn/winter (October- 
March, ONDJFM) and spring/summer (April-September, 
AMJJAS). There is precedent in the drought literature for 
an ONDJFM/ AMJJAS seasonal division (Kushnir et al. 
2010; Seager et al. 2005b; Seager 2007), and we found that 
our results and conclusions were relatively insensitive to 
the exact definition of our cold and warm seasons. 

2.2 Model simulations 

The model output comes from a series of simulations 
conducted with the NASA Goddard Institute for Space 
Studies (NASA GISS) atmosphere general circulation 
model, ModelE (Cook et al. 2009). In the SST-ONLY 
experiments, the authors forced a five-member ensemble 
(each member starting from different initial conditions) for 
the period 1932-1939 using only SSTs observed during 
this period. In the SST + LAND experiments, the authors 
ran an additional five member ensemble for the same 
period, incorporating not only observed SST information, 
but also estimates of dust aerosols and land cover changes. 
The land cover changes are integrated into the model by 
converting crop areas over the Great Plains to bare soil. 
The dust aerosols are included by adding a dust source over 
the Great Plains and allowing the model to deflate the dust 
into the atmosphere, where it can interact radiatively. 
Additional details on the experiments, including how the 
boundary forcing was estimated, are contained in Cook 
et al. (2009). For more detailed discussions of GISS 
ModelE formulations and performance, readers are referred 
to the available literature (Hansen et al. 2007; Miller et al. 
2006; Schmidt et al. 2006). Comparing the model preci- 
pitation anomalies against available observations, Cook 
et al. (2009) concluded that inclusion of surface changes 
and dust aerosols improved the simulated intensity, loca- 
tion, and spatial extent of the drought, compared to simu- 
lations with SST forcing only. At the time, however, a 
comparison of the circulation as characterized by geopo- 
tential heights and vertical velocities was not possible 
because of the lack of available data. 

To determine to what extent the dynamics in the 
COMPO reanalysis represent a response to SST forcing (as 
opposed to either noise or forcing from other factors), we 


also use a separate series of GCM experiments forced by 
idealized patterns of SST forcing in the tropical Pacific and 
Atlantic basins (Schubert et al. 2009). These simulations 
are described in more detail in Sect. 3. 

2.3 Other data 

Monthly gridded precipitation data are from the Climate 
Research Unit (CRU) at the University of East Anglia 
(Mitchell and Jones 2005). The CRU data are statistically 
interpolated from monthly station observations to a regular 
terrestrial grid at half degree spatial resolution and monthly 
temporal resolution, covering the time period 1901-2002. 
Gridded SST anomalies are taken from the HadISSTl.1 
dataset (Rayner et al. 2003), the same dataset used as part 
of the COMPO reanalysis and also used to force the 
aforementioned GISS model simulations (Cook et al. 
2009). Monthly climate indices for the PDO, AMO, and 
ENSO (NINO 3.4) are taken from various sources (Enfield 
et al. 2001; Mantua et al. 1997; Trenberth et al. 1997); 
these indices have been standardized (zero mean, unit 
standard deviation) to facilitate comparisons. 

3 Results 

3.1 Precipitation and SST anomalies 

Composite SST anomalies for the two drought periods 
show that both droughts were associated with warmer than 
normal north Atlantic SSTs and a cooler than normal 
eastern tropical Pacific (Fig. 1). Pacific anomalies are 
cooler and much more widespread during the 1950s, while 
north Atlantic anomalies are similar between the two 
droughts. The central north Pacific is anomalously warm 
during the 1950s, indicative of a negative phase PDO 
during this latter drought. Differences in the level of 
forcing can be seen even more clearly when climate indices 
during the droughts are averaged across years (Table 1). 
For example, NINO 3.4 and PDO index anomalies during 
ONDJFM are over twice the magnitude during the 1948- 
1957 drought, compared to 1932-1939. AMO anomalies 
are comparable, although the AMJJAS anomalies are 
stronger during 1932-1939 than 1948-1957. Dividing the 
Dust Bowl drought into two periods (1932-1934 and 1935— 
1939) shows the transient nature of the tropical Pacific 
forcing. NINO 3.4 anomalies were much larger during the 
earlier part of the drought, comparable to anomalies during 
the 1950s drought, and all but disappeared during the latter 
part. This implies that La Nina may have been important 
for the onset of the Dust Bowl drought (Seager et al. 2008), 
but may not be able to fully explain the drought persistence 
or the mean drought pattern. 
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Fig. 1 Background SST anomalies (°C) for 1932-1939 and 1948-1957, relative to 1921-1929. Anomalies for boreal winter (October-March, 
ONDJFM) and boreal summer (April-September, AMJJAS) are shown 


Table 1 Average NINO 3.4, AMO, and PDO index anomalies for 
boreal winter (ONDJFM) and summer (AMJJAS) seasons during the 
two droughts 



NINO 3.4 

AMO 

PDO 

ONDJFM 

1932-1934 

-0.48 

0.73 

0.04 

1935-1939 

-0.07 

0.79 

0.83 

1932-1939 

-0.23 

0.77 

0.53 

1948-1957 

-0.47 

0.75 

-1.11 

AMJJAS 

1932-1934 

-0.34 

0.89 

0.05 

1935-1939 

-0.20 

1.03 

0.60 

1932-1939 

-0.25 

0.98 

0.40 

1948-1957 

-0.35 

0.57 

-0.77 


All indices are standardized anomalies (zero mean, unit variance). 
Base period for standardizing is 1900-2009 for PDO and AMO 
indices. Base period is 1900-2008 for NINO 3.4 index because raw 
data for 2009 was incomplete. For the 1932-1939 drought, anomalies 
are also divided into two periods for the earlier (1932-1934) and later 
(1935-1939) periods of the drought 

The expected teleconnections of the Pacific and Atlantic 
ocean basins to NA climate in the COMPO reanalysis are 
well resolved, based on correlations between the climate 


indices and selected dynamical fields for 1909-1958 
(Fig. 2). During ONDJFM, the NINO 3.4 index is posi- 
tively correlated with 200 hPa heights in the tropics, 
flanked by negative correlations in the midlatitudes, espe- 
cially over the Pacific and Americas (Fig. 2, top). This 
represents the zonally and hemispherically symmetric 
response of the atmosphere to ENSO forcing alluded to in 
the introduction; the inverse pattern, associated with La 
Nina events, is associated with drought over NA. Addi- 
tionally, positive correlations with pressure vertical 
velocity at 500 hPa (indicating anomalous subsidence 
during El Nino events) are seen in expected regions: the 
Amazon, southern Africa, and the maritime continent (not 
shown). For the north Atlantic, the dominant season of 
influence is boreal spring and summer (AMJJAS). The 
AMO index is negatively correlated with 1,000 hPa heights 
over the Caribbean and subtropical north Atlantic (Fig. 2, 
bottom), indicating a weakening of the Azores subtropical 
high when the north Atlantic is warmer than normal (Sutton 
and Hodson 2005). The large area of positive correla- 
tions over the Pacific sector may indicate some influence 
from the tropical Pacific, since Atlantic variability, espe- 
cially in the tropics, is not completely independent from 
ENSO (Enfield and Mayer 1997). There is also a long band 
of negative correlations with pressure vertical velocity at 
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850 hPa, indicating the northward migration of the Inter- 
tropical Convergence Zone in the Atlantic sector when the 
North Atlantic is warmer than normal (Yoo and Zeng 
2009) (not shown). Teleconnections with the PDO during 
ONDJFM are quite similar to the NINO 3.4 correlations 
(not shown), again implying that the PDO may be forced 
wholly or partially from the tropical Pacific. Because of the 
difficulty in effectively differentiating PDO from ENSO 
climate impacts in the available datasets, we focus the 
majority of our analyses and discussion on climate 
responses to direct ENSO and Atlantic forcing. 

Both droughts were predominately warm season 
droughts, with the largest precipitation deficits and areal 
extents during AMJJAS (Fig. 3). To determine how well 
these drought patterns match the idealized pattern that would 
be expected from a combined Pacific Cold/ Atlantic Warm 
(PC AW) SST pattern, we compare them against an idealized 
drought pattern calculated from the CRU data using pattern 
correlations, the product-moment coefficient of linear cor- 
relation between variables corresponding to the same loca- 
tion. Pattern correlations are commonly used in verification 


of forecast fields and can be interpreted identically as a 
Pearson correlation. The PCAW pattern is calculated by 
averaging precipitation anomalies for all years, excluding 
the 1930s and 1950s droughts, with cold Pacific (—0.5 
standard deviation in NINO 3.4) and warm Atlantic (+0.5 
standard deviation in AMO index) SST conditions (1911, 
1918, 1944, 1945, 1976, 1989, 1999, 2000, 2001), maxi- 
mizing the SST-forced drought signal. Pattern correlations 
are calculated over the region 20°N-50°N and 130°W- 
80°W. This area encompasses the primary regions with 
significant teleconnections to Pacific and Atlantic SSTs, 
including the Pacific Northwest, the southwest, Mexico, and 
the southeast. PCAW 9-year composite SST and precipita- 
tion anomalies are shown in Figs. 4 and 5. 

Comparing Figs. 5 and 7, it is seen that the 1950s 
drought resembles a classic combined Atlantic and Pacific 
SST forced drought pattern, with drying extending in a 
band across the southwestern United States, Mexico, and 
into the southeastern United States, consistent with the 
forcing and teleconnections previously shown. The ideal- 
ized PCAW drought pattern has a much higher pattern 


Fig. 2 Teleconnections 
(Pearson correlations) between 
climate indices (NINO 3.4 and 
AMO) and dynamical fields 
from the early twentieth century 
reanalysis 


ONDJFM NIN03.4 Correlation: 200 hPa Height, 1909-1958 
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Fig. 3 Precipitation anomalies (mm day x ) for 1932-1939 and 1948-1957. Anomalies for boreal winter (October-March, ONDJFM) and boreal 
summer (April-September, AMJJAS) are shown. Data taken from the CRU 2.1 gridded dataset (Mitchell and Jones 2005) 


correlation (r = 0.58) with the 1950s drought than the 
1930s (r = 0.24) (Table 2). During the 1930s, the drought 
is centered predominately over the central Great Plains and 
has a much larger, continuous extent. The 1930s also show 
little or no drying in Mexico or the southwestern United 
States. As noted before, general circulation model experi- 
ments using observed SSTs from these periods are typically 
able to reproduce the pattern and intensity of the 1950s 
drought, but reproduce the intensity and spatial extent of 
the 1930s only superficially (e.g., Cook et al. 2009; 
Schubert et al. 2004a, b; Seager et al. 2008). 

3.2 Model responses to idealized SST forcing 

In order to determine to what extent the circulation 
anomalies during the 1932-1939 and 1948-1957 droughts 
fit with expected responses to SST forcing, we first 
examine results from a series of experiments using the 
Community Climate Model (CCM) version 3. These simu- 
lations are from an inter-model comparison conducted by 
the U.S. Climate Variability and Predictability working 
group, investigating the role of SST forcing and land- 
atmosphere feedbacks on drought (Schubert et al. 2009). 


In these experiments, the model was forced for fifty years 
using constant cold Pacific and warm Atlantic SSTs, 
essentially idealized versions of the boundary conditions 
present during the COMPO droughts, and compared 
against a run forced by climatological SSTs. In cases where 
the dynamics in the COMPO reanalysis match the idealized 
experiments, we may have some confidence that the 
drought pattern was SST forced. In cases where the 
COMPO dynamics and dynamics in the idealized model 
runs diverge, we search for an alternative explanation. 

Height and pressure vertical velocity anomalies in the 
idealized model experiments are shown in Figs. 6 and 7, 
respectively. Positive height anomalies extend from the 
north Pacific and across southwest and southeast NA, 
essentially the inverse pattern from the NINO 3.4 corre- 
lation plot in Fig. 2. There is also a clear wave train 
emanating from the tropical Pacific during ONDJFM, with 
negative height anomalies over the subtropical Pacific and 
northwestern NA and positive anomalies over the north 
Pacific and Greenland. During AMJJAS, height anomalies 
at 850 and 500 hPa over the subtropical Atlantic are nega- 
tive, consistent with AMO correlations from Fig. 2. 
Pressure vertical velocity anomalies are stronger during 
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Fig. 4 SST anomalies (°C) from the PCAW composited drought 
pattern (1911, 1918, 1944, 1945, 1976, 1989, 1999, 2000, 2001) 



Fig. 5 Precipitation anomalies (mm day x ) from the PCAW com- 
posited drought pattern (1911, 1918, 1944, 1945, 1976, 1989, 1999, 
2000, 2001) 


Table 2 Pattern correlations between observed (CRU) and modeled 
(SST-ONLY, SST + LAND) drought patterns over the region 20°N- 
50°N and 130°W-80°W 



CRU (PCAW) 

SST-ONLY 

SST + LAND 

CRU (1932-1939) 

0.24 

0.21 

0.41 

CRU (1948-1957) 

0.58 

0.23 

0.26 

CRU (PCAW) 

N/A 

0.06 

0.07 


PCAW refers to composited drought patterns from years with cold 
Pacific and warm Atlantic SSTs (1911, 1918, 1944, 1945, 1976, 1989, 
1999, 2000, 2001) 


ONDJFM, with anomalous subsidence in the Pacific and 
southwest/southeast NA flanking a band of anomalous 
uplift; though diminished, the subsidence over southwest/ 
southeast NA continues into AMJJAS. 

3.3 Observed droughts, height anomalies 
(200, 500, 850 hPa) 

In the COMPO reanalysis during ONDJFM (Fig. 8), both 
droughts have positive height anomalies in the north 
Pacific, especially at 200 and 500 hPa. For 1948-1957, this 
is indicative of a wave train emanating from the tropical 
Pacific, with positive anomalies in the north Pacific flanked 
by negative height anomalies over northwest NA and the 
subtropical Pacific. The wave pattern for 1948-1957 is 
shifted southwest relative to the idealized anomalies in 
Fig. 6; positive heights are over northeast Canada instead 
of Greenland, and the negative anomalies in northwestern 
NA extend too far to the southeast. Still, the 1948-1957 
anomaly pattern is largely consistent with the idealized 
response. Circulation anomalies during 1932-1939 diverge 
much more strongly from the idealized case in Fig. 6; 
positive height anomalies extend across NA zonally from 
the North Pacific to the North Atlantic basin. Both droughts 
lack the expected high pressure over southeast NA during 
ONDJFM, although 1948-1957 does show positive heights 
over the west and southwest at 200 and 500 hPa. 

In the warm season, AMJJAS (Fig. 9), the most notable 
deviation from the expected response is a strong ridge of 
high pressure at 200 and 500 hPa centered over central NA 
during 1932-1939. This feature is remarkably consistent 
with an independent analysis of upper air observations 
during the 1930s (Bronnimann et al. 2009), a study that 
also identified an upper level ridge over NA during this 
drought. For 1932-1939, there is also a widespread anti- 
cyclonic anomaly over western NA at 850 hPa which, 
when combined with the negative height anomalies over 
the southeast, acts counter to the climatological flow from 
the Gulf of Mexico and into the continent (see Sect. 3.5). 
Both droughts show the negative height anomalies over the 
tropical/subtropical north Atlantic during AMJJAS; these 
reflect a typical response to warm Atlantic SSTs (Fig. 6) 
(Kushnir et al. 2010). 

3.4 Observed droughts, pressure vertical velocity 
(500 hPa) 

Pressure vertical velocity anomalies at 500 hPa for both 
droughts (Fig. 10) correspond closely to the precipitation 
anomalies (Fig. 3) from the CRU data. For 1948-1957, the 
anomalous subsidence is focused over southwest and 
southeast NA, a typical response to the combination of La 
Nina forcing and a warm subtropical Atlantic. The 
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ONDJFM PAC-Cold/ATL-Warm Height Anomaly: 200 hPa 


AMJJAS PAC-Cold/ATL-Warm Height Anomaly: 200 hPa 



160°E 160°W 120°W 80°W 40°W 

ONDJFM PAC-Cold/ATL-Warm Height Anomaly: 500 hPa 




160°E 160°W 120°W 80°W 

AMJJAS PAC-Cold/ATL-Warm Height Anomaly: £ 




Fig. 6 Height anomalies (m) from CCM model experiments forced with an idealized cold Pacific and warm Atlantic SSTs. Anomalies are 
relative to a 50-year run forced by climatological SSTs 


subsidence during the 1950s persists through both seasons, 
remaining in roughly the same location but intensifying 
during AMJJAS. For 1932-1939, the anomalous subsi- 
dence also closely overlaps the precipitation anomalies, but 
is centered farther north over the central Great Plains. The 
1930s anomalies also have a much more widespread and 
continuous extent, spreading out to cover much of the 
continental United States. Notably, subsidence anomalies 
in the southwest and southeast are weaker for the 1930s 
compared to the 1950s, supporting the notion that SST 
forcing only weakly contributed to the spatial pattern and 
intensity of the Dust Bowl drought. As an aside, we note 
the close association between the subsidence and precipi- 
tation anomalies helps to increase our confidence in the 
validity of the COMPO reanalysis, since these two prod- 
ucts use completely independent underlying datasets. 

3.5 Observed droughts, Great Plains low-level jet 

An important feature of the warm season climate over the 
Great Plains is the Great Plains low-level jet (GPLLJ), 
characterized by a maximum in low level wind speeds, 
with the jet core extending from the Gulf of Mexico into 
the central Great Plains (Higgins et al. 1997; Weaver and 


Nigam 2008). The GPLLJ is strongest during May-June, 
although it is active throughout the warm season, and is an 
important mechanism for moisture transport from the Gulf 
into the continental interior (Weaver and Nigam 2008). 
The GPLLJ appears to have the strongest impact on pre- 
cipitation during late summer, July-September (Weaver 
et al. 2009). Variations in the strength of the GPLLJ have 
been linked to severe droughts and floods in the Plains in 
previous studies (Arritt et al. 1997; Mo et al. 1997). 

From the COMPO reanalysis, we calculated a GPLLJ 
index for AMJJAS by averaging meridional wind speed 
below 900 hPa from 25°N-35°N and 102°W-97°W (after 
Weaver and Nigam 2008), with anomalies calculated rela- 
tive to 1921-1929 (our common baseline period). The 
reanalysis shows marked variability in the strength of the 
jet (Fig. 11), with extended periods of low wind speeds 
during the 1930s (1932-1936) and 1940s (1942-1945). 
This is in sharp contrast to the 1948-1957 drought, when 
the jet is close to average strength or weakly enhanced. The 
weakening of the jet for the 1930s drought and the impact 
of the Atlantic on moisture fluxes from the Gulf of Mexico 
into the Great Plains has been noted previously (Bronni- 
mann et al. 2009; Schubert et al. 2004b). Our GPLLJ index 
was not significantly correlated to either Atlantic or Pacific 
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ONDJFM PAC-Cold/ATL-Warm Omega Anomaly: 500 hPa 



AMJJAS PAC-Cold/ATL-Warm Omega Anomaly: 500 hPa 



Fig. 7 Pressure vertical velocity anomalies (Pa s 1 x 1,000) from 
CCM model experiments forced with an idealized cold Pacific and 
warm Atlantic SSTs. Anomalies are relative to a 50-year run forced 
by climatological SSTs 

SSTs, which is not especially surprising because of the 
complex and seasonally dependent nature of correlations 
between the GPLLJ and SSTs (Weaver et al. 2009). 
Instead, the variability in the jet in the COMPO reanalysis 
appears to be tightly linked to dynamics over NA, as cor- 
relations with the geopotential height fields show (Fig. 12, 
left and middle panels). The jet weakens when high pres- 
sure dominates at high levels over central NA and at low 
levels over the west, and is associated with reduced 
atmospheric moisture in the central Great Plains region 
(Fig. 12, right panel). The relationship between the GPLLJ 
and circulation in our analyses is largely consistent with 
other studies (Weaver et al. 2009). 

As noted previously (Fig. 9), the 1932-1939 drought is 
dominated during the summer season by an upper level 
ridge and low level anticyclonic anomaly over western NA, 
features absent during the 1948-1957 drought. This 
appears to be the dynamical explanation for the weakening 
of the jet during the 1930s. The importance of the GPLLJ 
during the 1932-1939 drought can also be seen in the 
specific humidity anomalies (not shown). During the 
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1932-1939 drought there are strong reductions in atmo- 
spheric moisture over the central Great Plains; for 1948- 
1957 moisture reductions are confined predominately to the 
southwest and southeast regions, where SST forcing 
dominates. 

3.6 Model comparison 

Precipitation anomalies from a series of model experiments 
using GISS ModelE are shown in Fig. 13. In the left panel 
are ensemble mean anomalies from a 1932-1939 simula- 
tion forced by observed SSTs, with default land surface 
characteristics and no dust aerosols (SST-ONLY). On the 
right are anomalies from an ensemble for the same period, 
forced by observed SSTs, land degradation in the form of 
crop areas converted to bare soil, and radiatively active 
dust aerosols originating from a Great Plains dust source 
(SST + LAND). In the SST-ONLY case, the model sim- 
ulates a drought that closely resembles a typical La Nina 
type drought, with moderate drying in the southwest and 
southeast. When land degradation factors are included 
(SST + LAND), the drought intensifies and expands 
northward into the Great Plains. Pattern correlations of the 
modeled drought patterns against the observed precipita- 
tion pattern during the 1930s (Table 2) shows a large 
improvement in the simulated drought pattern when land 
degradation effects are included (r = 0.21 for SST-ONLY; 
r = 0.41 for SST + LAND). While results from the GISS 
model experiments are not directly comparable to either the 
1950s drought or the idealized (PC AW) pattern (due to 
different SST boundary conditions), these comparisons can 
still be informative. For example, addition of land degra- 
dation to the model (SST + LAND) does not improve the 
pattern correlation between the modeled 1930s drought and 
the observed 1950s and PC AW droughts. This fits well 
with existing evidence that these drought patterns are 
predominately tied to SST forcing, rather than forcing or 
feedbacks from the land surface. 

The additional drying in SST + LAND comes primarily 
from the dust aerosols, which reduce net radiation at the 
surface and the top of the atmosphere, primarily through 
enhanced reflection of shortwave radiation. The net radia- 
tion reductions drive a feedback similar to the ‘Charney’ 
mechanism (Charney 1975), where reduced energy inputs 
to the atmospheric column are compensated by subsidence 
and adiabatic heating in order to restore thermodynamic 
equilibrium. The anomalous subsidence, in turn, reduces 
precipitation primarily by suppressing convection (Cook 
et al. 2009). Vegetation losses act as a secondary amplifier 
by limiting evapotranspiration and moisture supply to the 
atmosphere, and lead to warming of the surface throughout 
the Great Plains as the Bowen ratio shifts from latent to 
sensible heating. Even with the improved precipitation in 
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ONDJFM Height Anomaly: 200 hPa, 1932-1939 ONDJFM Height Anomaly: 200 hPa, 1948-1957 








Fig. 8 ONDJFM Height anomalies (m) at 200, 500, and 850 hPa for 1932-1939 and 1948-1957. Anomalies are calculated relative to mean 
conditions for 1921-1929 


SST + LAND, there is still too much drying in Mexico 
and the southwest, but this may reflect either deficiencies in 
the GISS model simulation of ENSO teleconnections or 
other factors not accounted for in the model simulations. At 
the time when Cook et al. (2009) was published, large 
scale dynamical datasets such as the COMPO reanalysis 
were unavailable for comparison against the model simu- 
lations, and validation efforts were limited to comparisons 
against temperature and precipitation. With the COMPO 
reanalysis now available, we can revisit these model 
experiments with the purpose of determining if the dust 
aerosols and vegetation reductions modify the model 
dynamics in directions consistent with the COMPO 
reanalysis. 

The circulation response to land degradation in the 
model is strongest during late summer and early autumn 
(August-October, ASO). In the left column of Fig. 14 are 
model anomalies for precipitation, 850 hPa geopotential 
height, and 500 hPa pressure vertical velocity for the 
SST + LAND experiments. These are relative to a 1921- 
1929 model control run with SST forcing only, and are 
analogous to the anomaly plots we have seen previously for 
the COMPO reanalysis. In the right column of Fig. 14 are 


differences between the two 1932-1939 simulations 
(SST + LAND minus SST-ONLY). These represent the 
added effect of the land degradation (devegetation and dust 
aerosols) beyond the SST forcing. As seen in the precipi- 
tation anomalies, land degradation shifts the drought 
northward in the model, such that the drying is now 
focused over the northern and central Great Plains (Fig. 14, 
top panels). At low levels, a high pressure develops over 
central NA (Fig. 14, middle panels). This is especially 
clear in the difference plot on the right (SST + LAND 
minus SST-ONLY), where the addition of dust aerosols 
and land cover changes enhances the anticyclonic anomaly 
over western NA at 850 hPa. The anticyclonic anomaly 
acts counter to the climatological flow, weakening the 
model GPLLJ slightly in the SST + LAND experiment. 
For the SST-ONLY experiments, the 1932-1939 GPLLJ 
anomaly (relative to the model mean for 1921-1929) is 
slightly positive (+0.145 m s -1 ); in SST + LAND, the 
model simulates a slight weakening of the jet 
(—0.152 m s" 1 ). These anomalies represent departures of 
about half an interannual standard deviation for the 1921— 
1929 base period (0.364 m s -1 ), but are relatively small 
compared to the observed anomalies (Fig. 11). In general, 
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AMJJAS Height Anomaly: 200 hPa, 1932-1939 AMJJAS Height Anomaly: 200 hPa, 1948-1957 








Fig. 9 AMJJAS Height anomalies (m) at 200, 500, and 850 hPa for 1932-1939 and 1948-1957. Anomalies are calculated relative to mean 
conditions for 1921-1929 


the GISS model adequately reproduces the climatological 
magnitude of the GPLLJ, but has much diminished inter- 
annual variability compared to observations. 

The most important response in the model is the inten- 
sification and spreading of the anomalous subsidence 
(Fig. 14, bottom panels). The subsidence anomalies match 
well with the precipitation anomalies in the model, and are 
the primary driver behind the drying (via the ‘Charney’ 
mechanism, mentioned previously). The fact that the model 
can largely reproduce the drought without major changes in 
the GPLLJ suggests that the jet may be a secondary actor 
and may be most important for reinforcing precipitation 
anomalies that are controlled primarily by vertical motions. 
Dust aerosol effects in the GISS model are confined pri- 
marily to the lower portion of the atmosphere, below 500 
hPa, and do not act to enhance the upper level ridging 
observed in the COMPO reanalysis. 

4 Discussion and conclusions 

Dynamical studies of persistent drought events over NA 
(and the globe) have been limited by the availability of data 
(e.g., winds, pressure fields, etc). Over the last 60 years, for 


which comprehensive circulation data are readily available, 
there have only been two persistent North American 
droughts, 1948-1957 and 1998-2004. This has forced the 
drought research community to rely primarily on GCM 
experiments forced by observed SSTs to examine the 
underlying dynamics and to test our understanding of 
persistent NA hydroclimatic variability, especially for 
droughts prior to the 1948-1957 event (e.g., Herweijer 
et al. 2006, 2007; Nigam and Ruiz-Barradas 2006; Seager 
et al. 2005b; Seager 2007). While these modeling studies 
have been largely successful for almost all NA droughts, 
many models have difficulty simulating the Dust Bowl 
temperature and precipitation anomalies using SST forcing 
alone. This is not especially surprising, given that the Dust 
Bowl anomalies diverge widely from the typical pattern 
observed during other instrumental period droughts. Whe- 
ther this discrepancy with observations is the result of 
missing physical mechanisms (such as land degradation 
considered here), the superposition of unforced variance on 
top of the response to SST anomalies, or a consequence of 
model deficiencies has been difficult to investigate given 
the limited data available during the Dust Bowl. To help fill 
this research gap, we have used a new reanalysis product, 
combined with recent GCM experiments, to compare 
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ONDJFM Omega Anomaly: 500 hPa, 1932-1939 ONDJFM Omega Anomaly: 500 hPa, 1948-1957 



AMJJAS Omega Anomaly: 500 hPa, 1932-1939 AMJJAS Omega Anomaly: 500 hPa, 1948-1957 



Fig. 10 ONDJFM and AMJJAS pressure vertical velocity anomalies (Pa s 1 x 1,000) at 500 hPa for 1932-1939 and 1948-1957. Anomalies 
are calculated relative to mean conditions for 1921-1929 


AMJJAS Great Plains Low Level Jet Anomaly 



Year 


Fig. 11 Great Plains low-level jet index (ms x ) calculated from the 
COMPO reanalysis, averaged over AMJJAS. Index values are 
expressed as anomalies, relative to the mean for 1921-1929 


circulation features during the 1930s Dust Bowl and the 
1950s drought and offer support for a plausible mechanism 
to explain the differences: 

• To varying extents, circulation anomalies (especially in 
the tropical/subtropical Atlantic and north Pacific ocean 
basins) during the 1932-1939 and 1948-1957 droughts 
are consistent with forcing from a combined cold 
eastern tropical Pacific/warm Atlantic SST pattern. 

• Observed subsidence anomalies match well with pre- 
cipitation anomalies in both droughts, a noteworthy 
correspondence given that the two underlying data 


AMJJAS GPLLJ Correlation: 500 hPa Height, 1908-1958 


AMJJAS GPLLJ Correlation: 850 hPa Height, 1908-1958 


AMJJAS GPLLJ Correlation: 900 hPa Spec Hum, 1908-1958 
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Fig. 12 Pearson correlation between Great Plains low-level jet index (Fig. 11) and 500 hPa height, 850 hPa height, and 900 hPa specific 
humidity for AMJJAS 
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Fig. 13 Ensemble mean AMJJAS precipitation anomalies (mm SST + LAND includes observed SSTs and land degradation effects 
day -1 ) for two Dust Bowl drought simulations using GISS ModelE. (crop failure and dust aerosols) 

SST-ONLY refers to simulations forced by observed SSTs; 


products are based on independent observations. The 
precipitation and subsidence anomalies during the 
1950s occur predominately over the southwest, Mexico, 
and the southeast. For the 1930s, both precipitation and 
subsidence anomalies are centered over the central 
Great Plains, inconsistent with what would be expected 
from SST forcing alone. The idealized SST forced 
drought pattern (PCAW) correlates much more strongly 
with the 1950s drought (r = 0.58) than with the 1930s 
drought (r = 0.24). 

• The 1930s drought is dominated by an upper level 
(200 hPa) ridge and low level (850 hPa) anticyclonic 
anomaly over western NA, circulation features that 
appear to drive a sharp reduction in the strength of the 
GPLLJ. These patterns, and the weakening of the jet, 
are absent during 1948-1957. 

• Model simulations with dust aerosols and land degra- 
dation factors, in addition to SSTs, are able to improve 
the model simulation of some of the circulation 
anomalies, including the enhanced subsidence over 
the Great Plains and the low-level anticyclonic anomaly 
over western NA. However, the model is unable to 
reproduce the 200 hPa ridging and simulates a rela- 
tively minor weakening of the GPLLJ, even when the 
land degradation factors are included. 

Remaining differences between our simulations and the 
reanalysis could be related to either uncertainties in our 
estimation of the boundary forcing (e.g., magnitude of dust 
aerosol loading), other physical processes not represented 
in the model, or deficiencies in the model itself. Alternative 
explanations to explain the anomalous nature of the Dust 
Bowl drought have also been proffered in the literature. It 
has been speculated that the warm Atlantic may have been 
enough to extend the drought beyond the initial ENSO 
forcing, and that any differences between the Dust Bowl 


and other droughts, such as the 1948-1957 drought, can be 
attributed to internal, unforced variability (Hoerling et al. 
2009; Schubert et al. 2004a, b). Regarding the Atlantic 
influence, most model simulations of the Dust Bowl 
drought do include forcing from observed Atlantic SSTs 
(Hoerling et al. 2009; Schubert et al. 2004a, b; Seager 
et al. 2008), and recent work suggests that circulation 
responses to the Atlantic forcing are accurately represented 
in the models (Kushnir et al. 2010; Schubert et al. 2009). 
Even with presumably accurate representations of Atlantic 
and Pacific SST forcing, however, most models are only 
able to partially reproduce the pattern, intensity, and 
location of the Dust Bowl precipitation anomalies (Hoer- 
ling et al. 2009; Schubert et al. 2004a, b; Seager et al. 

2008) . In some multi-member ensemble simulations, pre- 
cipitation anomalies in certain ensemble members do bear 
a remarkable resemblance to the observed precipitation 
anomalies (Schubert et al. 2004a, b; Seager et al. 2008), 
supporting the claim that natural variability played a role. 
To our knowledge, however, even these SST forced simu- 
lations are unable to reproduce the continental scale 
warming associated with the Dust Bowl drought. This 
warming can be reproduced in model simulations when 
land surface degradation factors are included (Cook et al. 

2009) . Additionally, if natural variability was an important 
factor, one would expect to find Dust Bowl analogues in 
the extensive, 2000 year long paleo-record of drought for 
North America (Cook et al. 2004). Past droughts that bear 
a superficial resemblance to the Dust Bowl can be be 
identified; they are quite rare, however, and none fully 
capture the combined magnitude, duration, and geographi- 
cal extent of the Dust Bowl (Fye et al. 2003; Seager et al. 
2008). Seager et al. (2008), for example, found only three 
potential analogues, and even in those cases the pattern 
correlation match between those droughts and the Dust 
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ASO Precipitation Anomaly (SST+LAND) 
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Fig. 14 Anomalies (SST + LAND relative to a 1921-1929 SST 
forced run) and differences (SST + LAND minus SST-ONLY, 1932- 
1939) in 850 hPa heights (m), 500 hPa pressure vertical velocity 


ASO Prec Diff (SST+LAND minus SST-ONLY, 1932-1939) 
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ASO 850 hPa Height Diff (SST+LAND minus SST-ONLY) 



ASO 500 hPa Omega Diff (SST+LAND minus SST-ONLY) 



(Pa s 1 x 1,000), and precipitation (mm day x ) from GISS model 
simulations (SST + LAND minus SST-ONLY) 


Bowl was only 0.40-0.50. This suggests that natural vari- 
ability is insufficient, by itself, to explain the Dust Bowl 
drought. 

Natural variability and the land degradation effects may 
have both significantly contributed to the unique nature of 
the Dust Bowl, and this has been suggested by others 
(Bronnimann et al. 2009; Hoerling et al. 2009). The 
question remains open, but results and conclusions from 
the modeling studies with dust aerosols and land cover 
changes (Cook et al. 2008, 2009) are consistent with the 


COMPO reanalysis. Clarification of the importance of land 
degradation relative to natural variability during the Dust 
Bowl will require further model experiments, including 
better estimation of the dust aerosol loading and surface 
changes and inclusion of these effects in other GCMs. 
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